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Abstract 

This article reviews and extends recent results concerning entangle- 
ment and frustration in multipartite systems which have some symme- 
try with respect to the ordering of the particles. Starting point of the 
discussion are Bell inequalities: their relation to frustration in classical 
systems and their satisfaction for quantum states which have a symmetric 
extension. It is then discussed how more general global symmetries of 
multipartite systems constrain the entanglement between two neighbor- 
ing particles. We prove that maximal entanglement (measured in terms of 
the entanglement of formation) is always attained for the ground state of 
a certain nearest neighbor interaction Hamiltonian having the considered 
symmetry with the achievable amount of entanglement being a function 
of the ground state energy. Systems of Gaussian states, i.e. quantum 
harmonic oscillators, are investigated in more detail and the results are 
compared to what is known about ordered qubit systems. 



1 Introduction 

Entanglement is the type of correlations which can be shared only with a finite 
number of parties — or to express it with the words of C.H. Bennett: "Entangle- 
ment is monogamous" (cf. ^). This is maybe one of the main characteristics 
of entanglement and it clearly distinguishes entanglement from classical corre- 
lations. The present paper is devoted to investigate this "monogamy property" 
of entanglement in symmetric multipartite systems with the particular focus on 
the relation to frustration, the existence of local hidden variable models and 
to ground states of Hamiltonians having the considered symmetry. The arti- 
cle is based on a talk given at the QIT-EQIS workshop in Kyoto 2003 and it 
essentially reviews and extends results from |5] and 13 ^ |S] . 

In order to understand a characteristic feature of entanglement it is useful 
to study the counterpart in classical correlations first — this will clarify the dif- 
ference between the quantum and the classical world. To this end Sec|21will as 
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a starting point discuss frustration in classical systems and their relation to Bell 
inequalities. In Sec|2|we will then utilize these observations in order to prove in 
a simple way that there exists a local hidden variable model for quantum states 
having symmetric extensions. Sec01 considers more general symmetries and in- 
vestigates the relation between the problem of maximizing the entanglement 
between nearest neighbors under a global symmetry constraint (e.g. transla- 
tional symmetry) and the task of calculating ground state energies of certain 
nearest neighbor interaction Hamiltonians. Finally, Sec[Slwill apply the ideas 
of the preceding section to Gaussian systems with symmetries characterized by 
symmetric graphs. 

2 Bell inequalities and frustration in classical 
systems 

Constraints on the possible range of correlations in the form of inequalities have 
been investigated for many years, even before physicists developed an interest 
in that subject due to the work of Bell (H] (see the monograph by Frechet |7]). 
The relation between Bell inequalities and frustration in classical systems was 
then pointed out and studied in the early eighties in particular by Fine 

In spite of the simplicity of this connection there are, however, surprisingly 
many publications in the field of quantum information theory in which this 
knowledge has apparently disappeared. The following section will recall the 
relations between joint distributions, local hidden variable models, frustration 
and Bell inequalities. 

2.1 Frustration in classical systems is due to loops 

Let us consider a set of observables Xi, . . . , Xn described within classical prob- 
ability theory and let us denote the probability that observable Xi leads to the 
outcome x by Pxi{x). Obviously, there exists always a joint probability distribu- 
tion Pxi,...,x„ — Yii which returns all the single distributions as marginals. 
However, if we fix in addition the pair distributions Pxi,Xj for a certain subset 
of pairs (i, j) in a non-trivial way, a joint distribution with these marginals in 
general no longer exists. A necessary condition for the existence of is 
of course the compatibility of overlapping distributions in the sense that 

^Pxi,X:,(,xi,X2) = ^Pxi,x3(a;i,a;3) = Pxi(a;i). (1) 

However, it is in general not sufficient, and the standard counter-example is 
a set of three observables with three pair distributions corresponding to total 
anti-correlations. 

It is very useful to depict the problem graphically and to assign a vertex 
to every observable and an edge to each of the fixed pair distributions. In this 
picture, compatibility of the pair distributions is sufficient for the existence of 
a joint distribution Pxi,...,x„ if the graph does not contain any loops: 
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Figure 1: Frustration in classical systems is due to loops. Let the vertices 
and edges of a graph be assigned to observables and given pair distributions 
(correlations) respectivley. Then (a) for tree graphs there exists always a joint 
distribution returning all the given distributions as marginals, (b) loops may 
cause frustration and (c) a necessary and sufficient condition for the existence 
of a joint distribution in a two-colorable graph is given by Bell's inequalities. 



Proposition 1 Consider a graph, such that every vertex k corresponds to one of 
n observables Xk and every edge (k, I) corresponds to a given pair distribution 
PXk,Xi- If all overlapping distributions are compatible in the sense of Eq.^\) 
and the graph does not contain any loop, then there exists a joint probability 
distribution Pxi,....x„ which returns all the given distributions as marginals. 

Proof: The proof can easily be formulated in terms of an induction. Assume 
there exists already a joint distribution Pxi,...,Xi which is compatible with all 
fixed pair distributions Px^^Xi where k,l < i and has the same marginal Pxi as 
all given pair distributions Pxi.Xj where j > i. Then we can construct a new 
distribution 

Pxi,...,x^+i = Pxi,...,Xi Px^,Xi+i Px, (2) 

which has all the required properties for a set of ? + 1 vertices. Since each step in 
the induction adds one edge and one vertex to the graph corresponding to the 
previous step, the construction only works for tree graphs, i.e., graphs without 
loops — otherwise one would at some point have to add an edge without adding 
a new vertex. 

The final joint distribution is then given by 
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where ej is the number of edges at the j'th vertex and the products have to be 
taken over all edges {k,l) and vertices j respectively^. □ 

Prop, n] shows that classically frustration is due to loops. This is in contrast 
to quantum mechanical systems, which can be frustrated even without loops. 

^Note that the extension is not unique. Consider for instance a three- vertex graph 
with edges (1,2), (1,3) and joint probabiUty Pxi,X2,X3- In this case we could as well 



take Pxi,X2,X3 = 

Ea;3 Q{X1,X2,X3) ■■ 



- Xl ,X2 ;X^ 



■■ and P 



Q where Q is any function fulfilling 5Dx2 QC^li^ai^s) 
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+ Q>0. 
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2.2 Joint distributions and local hidden variable models 



If the graph characterizing the fixed correlations between pairs of observables is 
two- colorable^ hke in Fig^c) then a necessary and sufficient condition for the 
existence of a joint distribution is given by the complete set of Bell inequalities, 
i.e., by the existence of a local hidden variable model. In this sense the corre- 
lations of a quantum state which violates a Bell inequality lead to unresolvable 
frustration when described within classical probability theory. 

Let us briefly recall the above mentioned result. We say that a set of mea- 
sured correlations {PAtBj } between observables Ai, . . . , Am^ and Si, ... , 
admit a description within a local hidden variable model PIIHIEIj if we can write 

Pa^b, (a, b)= f M{d\)xA^ {a, A)xb, {b, A). (4) 

J A 

Here A £ A is the hidden variable and the source of the correlation experiment 
is characterized by the probabilities with which the different A occur, i.e., by a 
probability measure M on A. The response function XAiia, A) gives the prob- 
ability that measuring the system in state A with observable Ai leads to the 
outcome a. The locality assumption in Eq.l^J is expressed in the fact that the 
response functions factorize, such that XAi (a, A) does not depend on Bj and 
XBj {b, A) is independent of Ai . 

The local hidden variable description of a set of correlations is equivalent to 
the existence of a joint probability distribution 

Proposition 2 There exists a joint probability distribution for all given pair 
distributions {PaiBj} if and only if the correlations admit a description within 
a local hidden variable model. 

Proof: Let A ~ {Ai, . . . , v4„i^), a — (ai, . . . , a„j^) be vectors of observables 
and their respective outcomes, and similarly for B, b. If Pa.b is the joint 
distribution for all pair distributions Pa^.b^, then 

Pa„b, (a, /3) = X! ^A.Bia, b) Sa,ai Sp.bj (5) 

a, 6 

is an admissible (deterministic) local hidden variable model with Pa.b playing 
the role of the measure M and the two delta functions corresponding to the 
characteristic functions in Eq.Q). 

Conversely, if (A, M, {xA^}, {xSj}) are probability space, measure and 
response functions for a local hidden variable model, then 

PA.B{a,b)= I M{d\)l[xAAa^,X)ll^BAbJA) (6) 

«=i i=i 

graph is called two- colorable, bicolorable or bipartite if we can divide the set of vertices 
into two disjoint sets such that no two vertices within the same set are connected by an edge. 
This is equivalent to saying that all the cycles are of even length. 
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Figure 2: Consider a multipartite quantum systems composed out of a single 
particle (on the right) and (n—l) "layers" of m particles each. Assume complete 
permutation symmetry within each layer and let p„ be the n-partite reduced 
state containing one particle of each layer plus the remaining single particle. 
Then pn admits a local hidden variable description for all correlations of an 
arbitrary number of observables on the n'th site and m observables on each of 
the first (n — 1) sites. 



is a joint distribution, which returns all pair distributions PAi,B as marginals. 

□ 

Note that Prop[21 naturally generalizes to n-partite systems, where correla- 
tions of the form PAiBj...Nk instead of pair distributions are given. 

3 Quantum states with symmetric extensions 

As already mentioned frustration in quantum systems are possible even without 
loops in the configuration. The simplest example is a system of three qubits dis- 
tributed to Alice, Bob and Charlie, where Alice wants to be maximally entangled 
with both her colleagues. Clearly, this is not possible and one way to see this is 
to realize that in such a scenario Alice could teleport an unknown qubit perfectly 
to both. Bob and Charlie, which contradicts the no-cloning theorem fTD). In fact, 
if Alice is maximally entangled with Bob she cannot share any correlations with 
Charlie, neither entanglement nor classical correlations and it has recently been 
proven that there is even a quantitative trade-off between the amount of entan- 
glement Alice shares with Bob and the maximal amount of classical correlations 
shared with Charlie 

We will in the following restrict to symmetric situations, where we have for 
instance Alice being surrounded by several Bobs and Charlies in a symmetric 
manner. Using the observations of the previous section, we can easily recover 
results of jl| and 

Proposition 3 Let p be a density matrix acting on a Hilbert space TIa ® 
such that all bipartite reduced states of the form pABj are equal. Then the state 
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pAB = PABj admits a local hidden variable model with respect to all correla- 
tions between m observables on site B and an arbitrary number of observables 
measured on site A. 

Proof: Let {Fi{ai)}, {Gj{bj)} be the POVM operators corresponding to the 
measurement devices Ai and Bj with respective measurement outcomes Oi and 
bj. Then 



Pa., 



,iai,bi 



, bm) = tr 



pF^{ai)®I^Gj{bj) 



(7) 



is an admissible joint probabiHty distribution compatible with all pair distribu- 
tions PA,,Bj{ai,bj) = tr [pABFi{ai) (S) Gj{bj)]. By applying PropUland PropEl 
we have then that there exists a local hidden variable model for all correlations 
of an arbitrary number of observables acting on A and m measurement devices 
for B. □ 

Prop, inishows that if the number m of equal neighbors in a quantum system 
is increased, then the bipartite correlations between the central particle and one 
of the neighbors become more and more classical. In fact, is was conjectured by 
Schumacher and proven by Werner (cf @J|^]) that if m tends to infinity, then 
PAB cannot contain any entanglement — this is the only possibility of having a 
starlike extension with arbitrary many parties. 

The following result generalizes PropOlto cases in which local hidden variable 
models for more than two parties can be constructed: 

Proposition 4 Let p acting on 7i®^ be a density operator characterizing a 
state of N — {n — l)m + 1 particles with permutation symmetry within each of 
(n — 1) disjoint sets of m particles (see Fig\^. Take the n-partite reduced state 
Pn which contains one particle of each set plus the remaining particle. Then 
Pn admits a local hidden variable description for all correlations of an arbitrary 
number of observables on the n'th "single particle" site and m observables on 
each of the remaining (n — 1) sites. 

Proof: Let {-Fj'-'(aj*-')} be the POVM corresponding to the j'th observable 
on site i with respective outcome a^''. Then 



tr 



k=l i=l 



(8) 



is an admissible joint distribution for one observable (j) on the n'th site and m 
observables on the remaining (n—1) sites. Having the tensor factors properly 
ordered, this distribution is compatible with all correlations measured on p„. 
This is seen by noting that taking the sum over a[,*^ in Eq.® corresponds to 
tracing out the fc'th particle in the i'th set ("layer" in Figl^J. Due to the assumed 
permutation symmetry of p it does, however, not matter which (ni— 1) particles 
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in each layer are traced out — we will always end up with the n-partite reduced 
state pn- 

By PropHthe probability distribution in Eq.l(Sl) can be extended to a joint 
distribution including an arbitrary number of observables on the n'th site and 
by the multipartite generalization of PropI21 there exists a local hidden variable 
model for all the considered correlations.'^ □ 



4 Maximizing entanglement and minimizing en- 



Prop. 121 and Prop. 0] show that an increasing symmetry of a quantum state 
(in the sense of an increasing number of equal neighbors) constrains the entan- 
glement in such a way that violating Bell inequalities of a certain type is no 
longer possible. In the following we will consider multipartite quantum systems 
which have more general symmetries with respect to the ordering of the par- 
ticles and investigate quantitatively how this global symmetry constraints the 
entanglement between two neighboring particles. 

The considered symmetry group will always be a subgroup G of the group 
<S„ of all permutations of n parties. The global density operator p then com- 
mutes with all group elements^ and it is in particular invariant under the group 
averaging 



It is in many cases again advantageous to depict the problem graphically and to 
assume that the group G is the symmetry group of a graph whose vertices cor- 
respond to particles, i.e., tensor factors of the total Hilbert space. For simplicity 
we will only consider edge transitive graphs herc^. This provides a natural no- 
tion of "neighboring" particles, namely those corresponding to adjacent vertices 
in the graph, and we will not have to specify which neighbors we are consid- 
ering since they are all equal. Prominent examples of edge transitive graphs 
are stars, rings, cubic, hexagonal and trigonal lattices, permutational invariant 
clusters and the platonic solids. 

Typically, states with these symmetries appear as ground states (or equi- 
librium states) of particles on a lattice with equal nearest neighbor interaction 
along the edges of the considered graph^. We will show that the state which has 
the largest nearest neighbor entanglement under such a symmetry constraint 

•^Note that the conditions of Propl^can be weakened in two directions: First, permutation 
symmetry within each layer is not required — it is sufficient that all possible n-partite reduced 
states are equal to p„. Second, as mentioned in Ref.|SI, positivity of p is only required on 
product operators. Moreover, the numbers of particles per layer may be different. 

■'The group G is represented by unitary operators {Ug} which permute the tensor factors 
of the total Hilbert space H®^. 

^That is, every edge can be mapped onto every other edge by an element of the symmetry 
group. 

^if there is no symmetry breaking 



ergy 
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is always the ground state of a certain nearest neighbor interaction Hamilto- 
nian and that, moreover, the maximal achievable amount of entanglement is a 
function of the ground state energy. 

This statement is rather obvious if we quantify the entanglement using a 
linear functional like the the overlap / = {^-\pab\^-) with the singlet state 
of two qubits. In fact, in this case the maximum is always achieved for 
the ground state of the antiferromagnetic Heisenberg Hamiltonian, since we can 
write ^ 

$-)($-! = - 1 - ax ® (Jx ~ (Jy ^ (Jy ~ CTz ® (Tz ■ (10) 

In this way we can for instance relate the ground state energy density of the 
infinite antiferromagnetic Heisenberg spin-i chain to the maximal achievable 
singlet fraction of a state with infinite translation symmetry, which is then 
given by f^ax = ln27 

However, an analogous statement is still true if we measure the entanglement 
in terms of a highly nonlinear hmctional, the entanglement of formation\14:\. 
In contrast to the singlet fraction the entanglement of formation is a proper 
entanglement measure and it is closely related to the amount of pure state en- 
tanglement needed to prepare a state by means of local operations and classical 
communication. It is defined as 

Ef{pab) = inf |^Ki?(M',)| - p| , (11) 

where i?(^) — 5(tr^[|vl')(vl/|]) is the pure state entanglement of ^E" given by the 
von Neumann entropy S of the reduced state. Due to the fact that Ea. l)ll|l is a 
so-called convex hull construction, we can prove the following: 

Proposition 5 Consider the set Vq of multipartite quantum states with (per- 
mutation) symmetry group G of an edge transitive graph. The maximal entan- 
glement EpIpab) between two neighboring particles of a state p S Dg 'is attained 
for the ground state projector of a nearest neighbor interaction Hamiltonian with 
interactions along the edges of the considered graph. The maximal achievable 
entanglement is then a function of the ground state energy. 

Proof: First note that by the concavity of the entropy we can take the 
infimum in Ea. Hll|l over all decompositions of pab into mixed states {<Ti} as well. 
Moreover, Ep is by construction the convex hull of the functional x t-^ <S'(tr^ [x]) , 
which can equivalently be expressed as the supremum over all afhne functions, 
which lie below it (cf.^H])- That is 

Ef{pab) = inf l^S^piS{tTA['yi\)\^^PiCi = p) (12) 
= sup |tr [pABh] IVct : tr [ah] < SiiiAW]) \ (13) 



^Bounds on the maximally entangled fraction under symmetry constraints will be studied 
in greater detail in 1131 . 
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(14) 
(15) 



where h{s) is a fictive parametrization of the set of "interactions" h defined in 
Eo. lfTIjl . 

By assumption pab is a bipartite reduced state of a global state p G Vg 
which is invariant under averaging over the group G. We can therefore write 

tr [pabHs)] = tr [p{h{s) (g) l)] = tr [T{p) {h{s) ® l)] (16) 
= tr [pT{h{s) ® l)] =: tr [pH{s)] , (17) 

where H{s) is a "Hamiltonian" with equal nearest neighbor interactions h(s) 
along the edges of the considered graph. Due to the symmetry of H(s) we can 
now drop the constraint p e in calculating the maximal achievable Ef{pab)- 

sup Ep{pab) — — inf inf tr [Ti/(s)] = — inf eo(s). (18) 

Here, eo is the ground state energy of the Hamiltonian H{s) and the normalized 
projector onto the ground state space of the extremal Hamiltonian has then 
both, the required symmetry and the maximal entanglement properties. □ 

There are two drawbacks concerning the application of Prop|5| First of all, 
the parametrization of h{s) is in general not known explicitly. In fact, there 
are only very few systems for which we are able to calculate Ef{pab)- Besides 
highly symmetric one-parameter families of states E) this is at present only 
feasible for systems of two qubits ^Sl and symmetric two-mode Gaussian states 

m 

The second problem is the calculation of the ground state energy eo(s) for 
large systems. Though some particular two-qubit interactions lead to exactly 
solvable models in one dimension, the set of these models is apparently not large 
enough in order to answer the question about the maximal possible Ef{pab) for 
an infinite qubit chain in a straight forward manner. Using results from Ref . |2U| 
we can write the concurrence of two qubits (which is in turn a monotone function 
of £;f CHI) as 

c{pab) - max |o, - ^Jnf^^ tr [pab{X ® X^Y)] | , (19) 

where F is the flip operator^. For a chain of qubits this leads to a somehow 
deformed XXZ -I- Z-field model which is, unfortunately, not exactly solvable. 

However, for Hamiltonians which are quadratic in bosonic operators, i.e., 
interactions with Gaussian ground states, both problems — calculating Ep (that 
is parameterizing h[s)) and calculating eo(s) — are feasible. 

*The flip operator interchanges the two tensor factor: F|(/)) (X) = IV") ® \4>) 
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Figure 3: Examples of symmetric graphs are the five platonic sohds (e.g. the 
dodecahedron) and the hexagonal and trigonal lattice. If every vertex corre- 
sponds to a single mode of a Gaussian state then the entanglement Ep between 
nearest neighbors is maximized for the ground state of the nearest neighbor 
Hamiltonian given in Ea. H21|) . Although the entanglement is finite in each case, 
the respective ground states are infinitely squeezed. 



5 Gaussian states on symmetric graphs 

Let us consider a bosonic system of n modes described by a set of canonical 
operators (Qi, . . . , Qn, Pi, . ■ ■ , Pn) ='■ R obeying the canonical commutation re- 
lations [Qk,Pi] — iSkil- States with Gaussian Wigner distribution — so-called 
Gaussian states — are completely characterized by their first and second mo- 
ments with respect to the canonical operators [21]. Physically they may de- 
scribe modes of the electromagnetic field, atomic ensembles interacting with 
such fields, the motional state of a collection of ions in a trap or the low energy 
(bosonic) excitations of many other systems. All the information about correla- 
tions and entanglement properties of these states is contained in the covariance 
matrix 



Tki = tr 



p{Rk-{Rk),Ri-{Ri)} 



(20) 



where (Rk) = tr [pRk] and {•,■}+ is the anti-commutator. 

It was proven in Ref.jS] that for the case of Gaussian states on symmetric 
graphs^, where every vertex corresponds to a single mode (one quantum har- 
monic oscillator), the Hamiltonian whose ground state maximizes Ep^pAg) is 
of the form 

H - Y.^Qk + Qif + {Pk-Pif (21) 
(fc,;) 

= + 2^ 2^''''"' ^ ^^^^ 

(k,i) 

= : (23) 



^Symmetric graphs are those which are edge transitive and vertex transitive. Examples are 
rings, cubic, hexagonal and trigonal lattices, permutational invariant clusters and the platonic 
solids (see Figl^. 
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Ef{pab) 




Figure 4: The maximal achievable entanglement Ep (ebits) between two par- 
ticles in an A^-partite permutational symmetric cluster. The lower curve corre- 
sponds to A^-mode Gaussian states, and the upper curve shows the result for an 
7V-qubit system. Both vanish asymptotically as ~ log2 N. 



where the first two sums run over all edges {k,l) of the graph and the Hamil- 
tonian matrix h = hq ® hp is block diagonal with the blocks hq and hp being 
easily determined from the adjacency matrix of the graph. The ground state en- 
gj-gyio q£ ^Yyq Hamiltonian and the covariance matrix of the respective ground 
state To are then given by 

eo = \WhQhp\\i , Fo = hph^^ ® hghpK (24) 

Although the entanglement is finite in all non-trivial cases^^ all the ground 
states are infinitely squeezed, which is mathematically expressed in the fact that 
Tq has zero eigenvalues. 

Figs. |5l and 01 show the (analytic) results for a translational invariant ring 
and a permutational invariant cluster of harmonic oscillators and compare them 
to what is known from the case of qubit systems^^. Surprisingly, the maximal 
value of Ep is of the same order of magnitude for both systems and as it is 
expected for a ring of qubits, the ring of harmonic oscillators shows an odd-even 
oscillation with respect to the number of particles. 

^"The relation between eo and the maximal achievable entanglement Ep is in this case 
different than suggested in the proof of Prop|5] However, Ep is still given by a monotone 
decreasing function of eg. 

^^If we apply the interaction in Eg. 1211 to only two particles, then the ground state will be 
the infinitely entangled original EPR state 1221 . 

^^The analytic result for the permutational invariant qubit cluster was derived in Ref . 1231 . 
The case of qubit rings was investigated in Ref. 1241 . where a lower bound on the maximal 
nearest neighbor entanglement was derived. 
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Ef[pab) 




Figure 5: Maximal nearest neighbor entanglement (ebits) in a ring of N har- 
monic oscillators (the dotted curves represent the envelopes). In the limit 
N —> oo this approaches 0.30 ebits, which is comparable to the 0.29 ebits Woot- 
ters found as a lower bound for the case of an infinite qubit chain. 



This together with other results from Ref. 2_ indicates three different ten- 
dencies for the maximal Ep: 

1. It decreases with the number of adjacent vertices. 

2. It decreases with the total number of vertices. 

3. It is suppressed in loops with an odd number of vertices, which give thus 
rise to additional frustration. 
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